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1. Introduction. 

In this paper we consider a (2+l)-dimensional Toda-like chain 

Au n = (Vu n ,Vu n ) ( + ), (1.1) 

\U n -U n+ i U n -U n ^iJ 

where u n = u n (x,y), A and V are the two-dimensional Laplacian and gradient. This 
equation appears in the literature in various contexts. For example, it is known to 
describe the Backlund transformations of a Heisenberg-like magnetics (see pQ) and 
Laplace transformations of hydro dynamic-type systems in Riemann invariants [2]. In 
a recent paper [3] this system was shown to describe 0(3) a-fields (three-component 
vectors cr n of unit length) coupled by a nearest-neighbour Heisenberg-like interaction 
such as, e.g., graphite-like magnetics when the spins inside one layer are governed by 
the Landau-Lifshitz theory with an effective Heisenberg interaction between adjacent 
layers: H = H LL + "H H where 7i LL = J2 n with 

£ n = J dxdy (Vcr n ,Vcr n ) (1.2) 

and n H = \ En E P = n ±i Un P with 

Un P = g 2 J dxdy ln(l + (cr n ,cr p )y (1.3) 

A remarkable feature of the model (II. ip . which we will write in terms of the complex 
variables z = x + iy and z = x — iy, 

ddu n = (du n ) (Bu n ) ( + ) (1.4) 

\u n - u n+1 u n - U n -tJ 

where d and 8 stand for d/dz and d/dz, is that it can be viewed as a connecting link 
between almost all Toda-like chains. In the following section, we present its relationships 
with the Wu-Geng chain (WGC) [4], the 'two-dimensional Volterra equation' (2DVE) 
[5j E] , the two-dimensional Toda lattice (2DTL) [7] and the relativistic Toda chain (RTC) 

esi si ma mj 

The main goal of this work is to obtain some explicit solutions for (II. 4p . We will 
not elaborate the inverse scattering transform or the algebro-geometric approach from 
scratch. Instead, we use the links between our equation and the 2DTL together with 
the results [12] that give us possibility of bilinearizing (jl.4p and reducing it in section 
[3] to a system that belongs to the Ablowitz-Ladik hierarchy (ALH) [13J. Starting from 
the structure of the already known solutions for the ALH, we derive in sections H] and 
[5] Toeplitz and dark-soliton solutions directly from some determinantal identities. In 
section El we use the Fay identity for the ^-functions to derive the quasiperiodic solutions. 
Finally, in section [7] the results of sections HH6] are used to obtain the Toeplitz, dark- 
soliton and periodic solutions for the 2DVE. 
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2. WGC, 2DVE, 2DTL and RTC. 

The main equations of this paper can be viewed as a straightforward generalization of 
the WGC to (2 + 1) dimensions. Indeed, from equations (2) of jl], 

11 11 , n , 

Unt = , v nt = (2.1) 

V n +1 v n U n U n -\ 

where u nt stands for du n /dt, one can obtain that functions w n defined by v n = w n — w n -i 
satisfy 

w nt = h C, C = constant (2.2) 

Un 

which leads, after setting C = 0, to 

2 W n+ i - 2w n + U> n _i 

Wntt = w nt -. (2.3) 

(w n+ i - w n ) [w n - W n -i) 

Clearly, this equation coincides with fll.4p after replacing d and B with d/dt. 
On the other hand, rewriting equations (ll.4p in terms of the variables 

idu n -idu n 

o n = (2.4) 



(2.5) 



u n+l ~ U n U n — U n -\ 

one arrives at the system 

ida n = a n (b n+1 - b n ) 
idb n = b n (a n _i - a n ) 
which is closely related to the 2DTL: the quantities /„ defined by 

f n = \na n b n (2.6) 

satisfy 

ddf n = e fn+1 - 2e fn + e }n -\ (2.7) 

System (12. 5p is known since the works of Leznov, Savel'ev and Smirnov [6] where the 
authors demonstrated that this system, which was named the 'two-dimensional Volterra 
equation', represents the Backlund transformations for the 2DTL and constructed its 
general solutions in the finite case using the results of [HI IT5] . 

Another model that we would like to discuss is the RTC [SJ HI [TUl [11] which can be 
presented Hamiltonian system 

idq n = dn/dp n , idp n = -dUjdq n (2.8) 

with 

U = ePn {e gn+1 ~ qn - 1) • (2.9) 

n 

Equations ([2l§ 

' idq n = e Pn (e 9 "+ 1 - g " - 1) 

(2.10) 
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rewritten in terms of the variables 

a n = e p™+9n+i-9n > b n = e~ Pn (2.11) 

are 

idlna n = a n+ i - a n -i - l/b n+1 + l/b n (2 12) 

idlnb n = a n -i-a n . 

One can easily note that the second of the above equations coincides with the second 
equation of ( 12. 5ft . In a similar manner, starting from the Hamiltonian system 

idq n = dU/dpn, idp n = -dH/dq n (2.13) 

with 

U = e ~ Pn (e 9 "" 9 - 1 - 1) (2.14) 

n 

that has the form 

idq n = e~ Pn (1 — e 9 ™" 9 ™- 1 ) 

(2.15) 

one arrives, after using substituting (12. lip , at the system 

(2.16) 



id In a n = b n+1 - b n 



idlnb n = a n -xb n -ib n - a n b n b n+1 

whose first equation is the first one from ( 12. 5p . The relationship between the two RTCs, 
( I2.12p and (I2.16p . and our system (I2.5P becomes transparent if one considers RTCs from 
the zero-curvature viewpoint based on the scattering problem [91 [101 [TT] 

a n ^n+l - b~ = A - ll>n-l) ■ ( 2 -17) 

System (12. 12|) is the compatibility condition of (I2.17P and 

idip n = a n (ip n +i - tpn) (2-18) 
while (I2.16P plays the same role for (I2.17P and 

idif> n = b n (ij) n - 1 -ip n ). (2.19) 
As to our system, ( I2.5|) . it ensures the consistency of ("12. 18[) taken together with ( I2.19p . 

idj) n = a n (V> n +i - ^n) ^ 2 2Q ^ 

idlp n = b n {i> n -l ~ V'n) ■ 

To summarize, equations that are the subject of this paper can be viewed as describing 
the commutativity of two relativistic Toda flows. It should be noted that this 
commutativity leads, again, to the 2DTL. Indeed, considering equations ( I2.10p combined 
with ( I2.15P one can obtain by direct calculations that functions q n satisfy 

ddq n = exp (q n+ i - q n ) - exp (q n - g n _i) (2.21) 

which is another form of the 2DTL. 

Finally, comparing ( 12 .4p with (I2.20p one can conclude that solutions for our system, 
Uni are nothing but solutions for the auxiliary linear problems for the 2DVE; in other 
words, the ii n -model is dual to (a n , 6 n )-model. 
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(3.1) 



3. Bilinearization. 

To bilinearize our equation, that can be written as 

dBu n = p n (Un+l ~ 2u n + U n -l) 
(du n ) (Bu n ) = p n [U n+ i - U n ) (u n - U n -i) 

we start with the 2DTL, assuming 

T n — l^n+l /o ri\ 

Pn = o 3 - 2 

where solution for 

dd\nr n = Tn ~ lTn+1 (3.3) 

T 

n 

or, in the bilinear form, 

DD r n -r n = 2r n -iT n+1 (3.4) 

where D and D are the Hirota operators corresponding to d and d: D a-b = (da)b—a(db), 
Da-b= (da)b-a(db). 

One can easily obtain a wide range of solutions for the first equation of (13. ip by 
differentiating the 2DTL tau-function r n with respect to any parameter, 
d 

u n = — lnr„. (3.5) 

Since r n can be considered as a solution for all equations of the 2DTL hierarchy, i.e. as 
a function of an infinite number of 'times', r n = r n (£1,^2, •■■^1,^2, ■••)> ^1 = z i ^1 = Z, the 
last formula can be generalized as follows: 



3=1 



\ ~ 1 dlnr n _ dhlT n \ 



.1 



with arbitrary constants Cj and Cj. However, this apparently most straightforward 
approach is rather hard to implement because the second equation of (13.11) . which plays 
the role of the restriction for the set Cj and Cj, is a bilinear system which we cannot 
solve. That is why we will use another way to deal with our equations. The key point, 
that will be proved below, is that u n can be composed of two solutions for the 2DTL 
equation: 



u n = — (3.7) 



T, 
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where 



I DDu n ■ co n = w„_ia; n+ i + Xoo n (3.8a) 



\ DD T n -r n = r n _ir n+ i + A r z n (3.86) 

and A is a unimportant constant that can be eliminated by multiplying u n and r n by 
exp(— Xzz) without modifying u n . Calculating the first derivatives of u n , 

du n = \d u n ■ T n , du n = \d uj n ■ r n (3.9) 

T T 

n n 
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1 u n 

ddu n = -^DD u n -T n £ DD r n ■ r n (3.10) 

T T 
n ' n 

one can rewrite (13. ip as 

{DD — 2 A) U) n ■ T n = T n _iU n+ i + T n+ iU n _i (3 11^ 

The next step is to split the second equation of the above system in two bilinear 
equations. Introducing new r-functions f n , u n , f n and u n by 

iD uj n ■ r n = f n u n (3.12a) 
-iD u n ■ r n = f n u n (3.126) 

and noting that the r.h.s. the second equation of (13.111) can be presented as X n X n _$ 
where 

X n = T n uj n+ g — r n+ su n , 8 = ±1 (3.13) 
it is possible to achieve our goal by setting 

U n = U n -s, f n = f n+ s (3-14) 
which leads to 

iDu n -r n = f n _ s u n (3.15) 

—iDu n -T n = f n u n -s- 

Now we have to close this system adding the equations describing the dependence of the 
new r-functions on the variables z and z. These equations should be (1) bilinear, (2) 
compatible and (3) imply as a consequence the first equation of (13.111) . The resulting 
system can be written as the union of the 'positive' part, 

= iD co n ■ r n - f n _ & Gj n 

= iDf n -r n + f n - S r n+ 5 (3.16) 

= iD U) n ■ Cj n _s + U n -5&n, 

= iD u n ■ r n + f n u n _s 
= iD f n _ 5 • T n + f n r n _ 5 (3.17) 
= iD co n ■ u n + uj n+ suj n ^s, 
combined with the restriction 

= T n U n+ s — T n+ sOJ n — f n u n . (3.18) 

One can verify that these bilinear equations indeed lead to the solution of our problem. 
As to their compatibility (that can be checked directly, but after rather long and tedious 
calculations), we would like to note that they are a part of the extended ALH [16J. 
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Actually, we have passed from r-functions of the 2DTL to ones of the ALH using the 
results of [12] . In other words, the above proceeding can be viewed as an explanation 
of the following fact, that can be verified by straightforward algebra: the compatible 
system (I3.16p - (l3.18p which belongs to the extended ALH provides solutions for the 
system (13. ip . The reference to the ALH not only resolves the question of compatibility, 
but also gives us possibility of using the already known solutions for the Ablowitz-Ladik 
equations to build ones for the problem we are dealing with. 

Definitions ( 12.4D and the above equations yield the following expressions for 
solutions of the 2DVE: 

= 1: On — — ; > ° n = (3.19a) 

= -1 : a n = , b n = (3.196) 

T n UJ n +l T n UJ n 

that will be used in section [7J 



4. Toeplitz solutions. 



This type of solutions can be constructed of the Toeplitz determinants 

A^ = det|a fe+0 _ 6 | a)6=w . (4.1) 

The main idea is that among various determinantal identities one can find ones that 
are similar to the equations we want to solve. In the appendix we derive the necessary 
formulae by applying the Jacobi identity to the determinants (14. ip and the framed ones, 

i C C 2 ■■■ c £ 

a k a k+ i . . . a k +t-i 



F fe 



t+i 



(0 



(4.2) 



(see (IA.28p - (1A.33p ). These formulae can then be rewritten in terms of the border 
determinants B|, similar to (14. ip . instead of the (£ + l)-order ones, F^ +1 (C), 

l(C) = (-O'b^CO, 



where 



with 



B, fe (C) = det|/3 fc+a _ fe (C)| a , b=1 , 



The set of identities that we need to solve our equations is 

= D + • A k e - C^ttl B "-i 1 

= D + B k e ■ A k e+1 + A k tl B^- 1 

A k+L R fe 
D i-i 



= D+ B» • A k+1 Afe+1 Dk 



(4.3) 

(4.4) 

(4.5) 

(4.6a) 
(4.66) 
(4.6c) 
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= D_ B{ ■ A{ + C^Ati B^j (4.7a) 

= D_ Bj ■ A, fc+1 + Aj +1 B^ 1 (4.76) 

— D_ ■ A^ 1 + fk\ +l B^ +1 (4.7c) 

together with 

= A, fc+1 B£ - A, fc B, fc+1 + C^ttl^t-v (4-8) 

Here, D± stand for the Hirota operators corresponding to d± and d±a.k = otk±i- 

Now, the problem of finding solutions for (I3.16p — (13. 18|) becomes a 'combinatorial' 
one: one has to select r n , u n , f n , u n , from A*, B^ with different k and £. Below, we 
obtain three families of solutions: infinite (with respect to n), semi-infinite and finite 
ones. 

4-1. Infinite chain. 

This type of solutions corresponds to the following choice of the r-functions: 

r n = r* A?, u n = Wl (4.9) 

and 

f n = % BJL X , U n = NlXl (4.10) 

with constant From (14.81) . one can immediately derive that 

= r n u n+1 - T n+l U n - f n u n (4.11) 

provided 

t*lu* = (T*a)*. (4-12) 

Thus, our r-functions solve (13.181) with 5 = 1. Furthermore, equations (I4.6al) - (l4.7cl) 
lead to 

= D + u n ■ r n - f n _i Gj n (4.13a) 

= D + f n ■ r n + f n _x r n+ i (4. 136) 

= -D+a;„ ■ a) n _i + a; n _i a) n (4.13c) 

and 

= £>_ w n • r n + w n _i f n (4.14a) 

= L>_ f n _! • r n + r n _! f n (4.146) 

= D_u) n ■ u n + u n -iu n+ x. (4.14c) 
It is clear that to complete solution for ( 13. 16[) — ( 13. 1 7[) one has to meet 

d + = id, d. = id (4.15) 
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or, to take a k to be solutions for the linear system 

m k = a k+1 
ida k = a fc _i. 

The 'symmetric' set of solutions can be obtained by using, instead of (I4.3p . another 
representation of the determinants F^ +1 (£): 

Fj + i(C) = Q fc (C) (4.17) 



where 



with 



Q fc (C) = det| 7fc+a _ b (C)| a , = , , (4.18) 



7fc(C) = "fc - C«fc-i- ( 4 - 19 ) 

The calculations similar to the ones presented above demonstrate that r-functions 
defined by 



and 



with 



T n = n A", 0J n = LU* Q (4.20) 

r„ = f»C2+ 1 , w n =w»A? +1 (4.21) 
C^w* = — f*u* (4.22) 



solve flXT6|) - (lXT8|) with 5=1. 

Both these sets of solutions can be written as 

det larL. h (z, z)| , 1 . 

/ -\ 1 n+a-o \ ' / \a,b=l,...,£ /, OQ \ 

u n {z,z)=u* — — — — (4.23) 

det|a n+a _ fe (z,z)L,&=w 

where £ is an arbitrary positive integer and the elements of the determinants are given 

by 

a k (z,z) = J dhce(h) h k exp[—iQ h (z,z)) (4.24a) 

af {z, z) = J dh a(h) [1 - (h/Cf 1 ] h k exp [-iG h (z, z)\ (4.245) 

(we replaced f3 k and ^ k with ct k ) with arbitrary contour T, function a(h) and constant 
it*. The 'dispersion law' 9^ (z, z) is given by 

S h (z,z) = hz + h~ l z. (4.25) 

As an example, let us consider one of the simplest solutions of (I4.16p . Noting that 
system (14.161) leads to the Helmholtz equation, dda k + a k = 0, and rewriting the latter 
using the polar coordinates, z = re t9 , one can obtain in a standard way 

a k = exp[-i{6 + n/2)k]J k {2r) (4.26) 
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where is the kth Bessel function. These solutions correspond to (I4.24al) and (14.2461) 
with a(h) = (27iih)~ 1 and T being the unit circumference: T = {h : \h\ = 1}. In the 
'elementary' case of I = 1 expression ( 14. 23ft can be written as 

u n = u* + u ± ^f^e^ (4.27) 

where and u ± (u± = im*^ 1 ) are arbitrary constants. It is easy to see that solutions 
we have obtained are complex and singular, which, however, does not mean that they 
are non-physical. Say in the case of [3J, that was discussed in the introduction, the spin 
components, cr n = (an\ cr n 2 \ cr n ) , are related to u n by an + i&rP = 2w n /(l + \u n \ 2 ) 
and = (1 — |w n | 2 )/(l + |w„| 2 ). Thus, in the general case (a^ ^ 0) u n is 
complex and singularities of u n correspond to the vertical (southward) orientation of 
a n : u n = oo ^ = a n 2 ^ = and = —1. 

4-2. Semi-infinite chain. 

This type of solutions appears if one identify the index n with the size of the determinants 
A*, B|. Consider the functions T n , W n , f n , W n given by 

T n = T*A k , W n = W*B k (4.28) 

and 

f n = t B k n ~\ W n = W* A k n X\ (4.29) 
with constants T*, W*, f*, W* being related by 

T m W. = -£f.W.. (4.30) 
It is straightforward to verify that they solve 

= T n W n+1 - T n+1 W n - f n W n (4.31) 



as well as 



= D + W n -T n + T n ^W n 

= D + f n -T n -f n „iT n+l (4.32) 
= D + W n ■ W n . x - W n _ x W n 



and 



= (£)_ + C" 1 ) W n • T n + W n - 1 f n 

= {D_ + (- 1 )f n _ 1 -T n + T n _ x f n (4.33) 
= (£)_ + r x ) W n ■ W n + W n ^ W n+1 
that almost coincide with H3. 16|) — A3. 18|) for 5 = 1 after identifying 

d = id + , 5 = -id-. (4.34) 

The extra terms can be eliminated by exp(iz/() factor and one arrives at the following 
solutions: 

r n = r* exp(^/C) A k , u n = uj* B k n (4.35) 
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and 



exp(iz/() A 



fe+i 

n+1 



where Tjf , cu* , , cu* are arbitrary constants and 



c 



(4.36) 



(4.37) 



In a similar way one can derive the 'symmetric' set of solutions using the 
determinants C*. These two sets of the Toeplitz solutions u n can be written as 

, s det \oti, n h (z, z) I , , 
u n (z, z) = - J k+a - b ) _ (4.38) 

det|a fe+a _ fe (z, z)\ a ^ b=1 ^ n 

with an arbitrary positive integer k. The elements of the determinants are given again 
by (I4.24a|) and (14.24 6|) with the 'dispersion law' 

h (z, z) = -hz + hr x z (4.39) 

while the phases <p± are given by 

<t> + {z,z) = -z/(, </>.{z 1 z) = Cz. (4.40) 

The above formulae for the Toeplitz solutions demonstrate the typical for integrable 
systems phenomenon: solutions for the nonlinear equations are determinants of matrices 
that satisfy linear ones. This rule in our case can be extended as follows: the 
relationships between the r-functions r n and u n (which are, recall, solutions for the 
2DTL used to construct u n ) become linear when rewritten in terms of the 'inside- 
determinant' objects, otk and af. 

As in the case of the infinite chain, let us consider one of the simplest solutions for 
(I4.34p . Using the polar coordinates z = re ld one can obtain 

a k = exp [-i(9 + ir/2)k] 4(2r) (4.41) 

where Ik are the modified Bessel functions. This leads, together with definition 1 14. 40 p 

of 0±, 



to 



Mr, 



± 



ajf = exp [-i(0 + 7r/2)fc] 
and consequently to 



4(2r) + 



Mr, 



i r 



<fc±i 



(2r) 



(4.42) 
(4.43) 



det 



icf>±(r,6)_ 



7 fc+a _ 6 (2r) + ^/ fe+0 _ 6±1 (27 



a,6=l,...,n 



(4.44) 

det|/ fc+a _ 6 (2r)| a>6=1 _ n 

In context of the Heisenberg-like model [3] these solutions describe circular magnetic 
domain structures. 
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4-3. Finite chain. 

The semi-infinite solutions derived in the previous subsection can be easily modified to 
provide the finite ones with u n ^ for n = 0, 1, N only. To this end one has to define 

u {z,z) = uj^ z >*\ (4.45) 

where <f>± are defined by (14.401) . and to replace the integrals in (I4.24q|) and (I4.246j) with 
finite sums, 

TV 

a k (z, z) = ^ & p K ex P [~ iQ h p {z, z)] (4.46a) 
P =i 

N 

of (z, z) = ®p I 1 - ( VCf 1 ] ^ exp [-iQ hp (z, z)\ (4.466) 
P =i 

where 

Q hp (z, z) = -h p z + h~ l z. (4.47) 

These modifications do not change the fact that ida k = — a^+i and ida k = ctk-i which 
implies that u and u n given by ( I4.38P still solve ( II. 4p for n = 1, ...,N — 1. Thus, it 
remains to prove that they solve (11.4j) for n = 0, iV as well. The case n = is trivial 
because for both choices of (fi± 

ddu = (du )(du ) = (4.48) 

converting equation (ll.4p with n = into trivial one. Considering the right-end 
equation, one can show, using the identity 



det 

that 



N 

a—b 



E 

P =i 



N 

a 



(-D^'n^T n iXp _ XtfA4Ag) 



a,b=l,...,N P= l Xp 1<P<1<N 



u N {zrz) = <e l ^ z ~ z \ (4.50) 

where u[ is a constant: the factors exp [— i®h p ( z , z)\ cancel themselves when we divide 
det |a^ +a _ b | by det | ct fc _ 1 _ a ._ 6 1 leading to 

N 

tt: = «.iJ[ 1 -(VC) ±1 ] (4.5i) 

P =i 

Hence, 

ddu N = (du N ){du N ) = 0, (4.52) 

ensuring solution of f ll .41) for n = N. 

To summarize, N + 1 functions, Uq given by (I4.45P and 

det \a^ +a _ b (z, z)\ ab=1 

u n {z,z) =u (z,z) — — — — — — , n = l,...,N (4.53) 

det \a k+a - b {z,z)\ ajb=lj jn 

with ( I4.46ap -( 14T4T|) solve the system of iV + 1 equations ( TT4) for n = 0, N. 
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The dark-soliton solutions for our equations can be constructed of the determinants 

ft (A) = det 1 1 + A | (5.1) 

of the N x N matrices A that satisfy the 'almost rank-1' condition 

LA - AR = \£}(a\. (5.2) 

Here 1 is the N x N unit matrix, L and R are the constant diagonal matrices, | £ ) is 
the constant iV-component column, \ £) — (£ x , ...,£^) T and (a\ is the A-component row 
depending on the coordinates, (a (z, z)\ = (a\ (z, z) , ajv (z, z)). The second part of 
the 'solitonic ansatz', except for f 1 5 . 2 j) . is that the dependence of A on all coordinates 
(z, z and n) can be built by means of the shifts 

T c ft = ft (AH f ) (5.3) 

with the matrices being defined by 

Hc = (L-C)(R-C) _1 (5.4) 

(we do not write the unit matrix explicitly, so (L — Q stands for (L — £1) etc). 
The remarkable property of the above matrices is that the determinants 

ft c = T c ft, ft €r , = T^T^ft (5.5) 

satisfy the Fay-like identity 

{i - rj) ft c Qtr, + (v-C) n v t + (C - a* % = (5.6) 

(see, e.g., appendix of p2j). Introducing the differential operators <9g defined by 

T^T^ft = ft + i8d^Q + (5 2 ) (5.7) 

one can derive from (15.61) the differential Fay identities 

(f - a) iDt. ft Q • ft = (T^ftJ (T € ft) - ft a ft (5.8) 

and 

(£ - <*)(£ " 0) iD e ft Q • ft^ = (a - (3) [(T^ft a/3 ) (T 5 ft) - ft^] . (5.9) 
Applying these identities two times, one arrives at 

X - (£ - r,) 2 D^D V ft • ft = ft 2 - (T^T^ft) (T-%ft) (5.10) 

from which it is easy to obtain solutions for the 2DTL by associating the d, 8 operators 
with <9g and d v (with fixed £ and rj) and introducing the n-variable by means of the 
powers of T v 

ftn = T^ft (5.11) 

with 

H 5 = H„H„. (5.12) 
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From the Fay identities (15. 6p and (15. 9p one can derive after straightforward calculations 
that the functions 

T n =T*h n a T n u Q a (5.13a) 

W n = W^T n u % (5.136) 

f n =f,h n T n u T- l ^ a p (5.13c) 

W n = W*h n T n v ^ (5.13 d) 
where the /i-factors are defined by 

K = , hp = -, (5.14) 

r] — a 7] — p 

h = fih a hp(£-Ti), h = p,(£-r)) (5.15) 

(with constants \i and p, satisfying fip, = (£ — r?)~ 2 ) and the constants T*, W*, T*, W* 
are related by 

T*W* = T*W* — r (5.16) 

\a-VKP-V) 

satisfy the following equations: 

[zD e + a)} W n ■ T n = ixt n _ x W n (5.17a) 

[iD^ + \e(P,7))]f n -T n = - ntn-xTn+x (5.176) 

[iD^ + X^, V )]W n -W n ^^ -iiW n . x W n (5.17c) 

and 

[iD v + X v {/3, a)} W n ■ T n =JiW n . l f n (5.18 a) 

[zD„ + \ v {(3, £)] f n _! • T n = liT n _ x f n (5.186) 

[i^ + V^OlW-n-Wn =/iW n -iW n+ i. (5.18c) 

A e (a,/3)=A € (a)-A € (/3) (5.19) 

Comparing the above equations with (13 . 16[) and (I3.17P it is easy to conclude that 
to obtain solutions for our problem one has to identify 

d = fT^t 3 = -/Z -1 ^ (5.21) 

and to introduce some linear in z and z phases to eliminate extra A-terms in (j5.17aj) - 
( 15. 18 cj) . This leads to the following expression for u n : 

u n {z, z) = h» (5.22) 

det 1 + A n (z, z) Hq, 



Here 



with 
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where 

h = ^, (5.23) 

h a 

A n (z, z) — A (z, z) H™ (5.24) 

and 

<p = fi^ 1 X^(a, (3)z + fx' 1 A r) (/3, a)z + constant. (5.25) 

Finally, it remains to resolve the restriction (I5.12p and to write down explicitly the 
dependence of A on the coordinates, which can be done by calculating the limit in (15. 7h . 

iA~ x d^A = (L - R) (L - Cr 1 (R - C) _1 (5.26) 

(for any (). 

The restriction (I5.12p implies that the matrices L and R are not independent, 

(L-0(R-0 = ("-0(l-Ol- (5-27) 
Introducing matrices L and R, 



satisfying 



L = — !— (L-i/), R = -^—(R-u) (5.28) 

L-l) (R-l) =|^1 (5.29) 
and parameter /, instead of \i and ft, 

f= ^_ j =m~ri) (5-30) 
the z- and ^-dependence of A and <fi can be presented in a symmetric form as 

A(z,z) = A*exp{iM(z,z)} (5.31) 
where A* is a constant matrix, 

M (z, z) = f (R - L\z + r 1 (R- 1 - IT 1 ) z, (5.32) 

and 

(z, z) = f (h/ - h- 1 ) z + f- 1 {h p - h a ) z. (5.33) 

These formulae together with (15.221) describe the A-dark-soliton solutions for our 
problem. 

The simplest (one-soliton) solution can be written as 

u n (z, z) = Uif e 2ilpn{z ~ z) [1 + tanhp tanh[/ n (z, z)} (5.34) 

where 

p = iln 1 -^ 1 -^, (5.35) 
2 1 — h a R 1 — hgL 
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(note that now L and R are just complex numbers) the phase ip n is given by 

<p n (z,z) = £ z + r] z + ( n (5.36) 

with 

£o = { {hp 1 - hZ 1 ) , r) = ^(hp- ha) (5.37) 

and 

sin 2 Co = £o?7o, (5.38) 
while U n is given by 

U n (z,z) = £ s z + r] s z + ( s n (5.39) 

with 

& = ^ (R - L) , Vs = ^j (iT 1 - L- 1 ) (5.40) 

and 

sinh 2 C s = 6^. (5.41) 

To conclude, we would like to note that in the soliton case the link between r n 
and u n is linear in the terms of the matrices A n : the matrices that appear in the 
determinants in (15.221) are related by the constant diagonal matrix H/jH" 1 . 

6. Quasiperiodic solutions. 

In this section, we derive the periodic solutions for our equation proceeding in the way 
similar to the one used in the previous section. The main difference is in the starting 
point: instead of identity (15. 6p we use the original Fay's trisecant identity (see (16.71) 
below) . 

The solutions that we derive below are combinations of the ^-functions defined over 
a compact Riemann surface T of the genus g for which one can choose in a standard 
way a set of closed contours (cycles) {a«, &i}i=i,..., g with the intersection indices 

a.; o aj = bi o bj =0, cii o bj = 5ij i, j = 1, . . . , g (6.1) 

and g independent holomorphic differentials Wk satisfying the normalization conditions 



w k = Sik, i,k = l,...,g. (6.2) 
The matrix of the 6-periods, 

rt ik =(pzu k , i,k = l,...,g (6.3) 



determines the so-called period lattice, Lq = {m + Qn, m,n G Z 9 }, and the Abel 
mapping from T to the 2g torus C 9 /Lq (the Jacobian of this surface), 

P-> I u (6.4) 

JPo 
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where P is a point of T, u) is the g-vector of the 1-forms, u = (wi, . . . , zu g ) T , and P is 
some fixed point of T. 

The ^-function, 9(C) = 6(C,ty, is defined by 

0(0 = Yl ex p{™ (n,ttn) +2ni(n,C\} ( 6 - 5 ) 



where (n, C) stands for Y2i=i n id- This is a quasiperiodic function on C 9 

6(C + n) =9(0 (6.6a) 
9 (C + fin) = exp (-iri (n, ttn) - 2m (n, c) } 9 (C) (6.66) 

for any neZ 9 . 

The calculations presented below are based on the famous Fay's trisecant formula 
[T8| [19] that can be written as 

4 el 9 P A 9% - e p A e% 9 P 9% + e p e B 9 9 P AE \ = 0. (6.7) 

Here 

9%::?: = 9 U +Y1 J p "j ( 6 - 8 ) 

and the skew-symmetric function Bp, Ep = —£q, is given by 

e Q P = 9(e + J u) (6.9) 

where e is a zero of the ^-function: 9 (e) = 0. 

Now let us define the differential operators d x by 

9 p x = 9 + e p x d x 9 + o(e p x ). (6.10) 

In what follows, we use Oa and Ob defined near two points A and B that are fixed. One 
of them will play the role of d and another of d. Taking the limit in (16. 7p one can obtain 
the differential Fay's identity 

[D x + X X (P, Q)\ 9 p -9 = lx (P, Q) 9 P 9* (6.11) 

where D x is the Hirota operator corresponding to dx and the functions \x(P) and 
lx(P, Q) are defined by 

\ X (P,Q) = \ X (P)-\ X (Q) (6.12) 

with 



X X (P) = lim £ \ f (6.13) 



and 



y^x e Y x e p x 



e p 



7x(P,Q) = ^7T. (6.14) 



£ v £ 



X °x 
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Consider the functions 

T n =T*[h{P)] n e(t n + (6.15a) 

W n = W. [h(Q)] n U n + £ <A (6-156) 

f n =th n e(c n + [ Q ") (6.15c) 



V JAB J 

W n = WJi n 9{Cn) (6.150 



where 



h{P) = 4 (6.16) 

£ B 

and the n-dependence is given by 

Cn+l = Cn + / (6.17) 

It follows from (16. 7p that if one imposes the restrictions 

£ P £ B 

tw = T*W* ( 6 - 18 ) 

£ B £ B 

and 

hh = h{P)h{Q), (6.19) 
then these functions satisfy 

T n W n+l - T n+l W n = f n W n . (6.20) 
Furthermore, by taking 

where \i A and /ig are two constants related by 

HAUB (ef) 2 = 1 (6.22) 

one can get from (16. lip the set of identities which will be associated with the <9-fiows: 
[D A + X A (Q,P)}W n -T n = -fx A f n ^W n (6.23a) 
[D A + X A {Q, B)\ f n ■ T n = ft A f n ^T n+1 (6.23b) 

[D A + X A (Q,B)}W n -W n -i = ii A W n _ x W n (6.23c) 
and another one, 

[D B + X B (Q, P)) W n -T n = -pL B W n -i f n (6.24a) 

[D B + X B {Q, A)} r n _! -T n = -fi B T n _! f n (6.246) 

[D B + X B (Q, A)} W n ■ W n = - fi B W n -x W n+1 (6.24c) 
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associated with the (^-equations. 

Comparing the above equations with (13.161) and f)3.17p . one arrives at 

i - 1 

d = — d A , d= - — d B 
Ha m B 

which determines the dependence on z and z, 

C n (z, z) = + za + zb + nc 
where c was defined above, 



B 



while 



a = if lim 



(jO 



i f- 1 lim 4 / 
p^ b £ B Jb 



where we have introduced the constant / instead of ha and hb, 



f - ——j ~ »b£ b - 



(6.25) 
(6.26) 
(6.27) 

(6.28a) 
(6.286) 

(6.29) 



Again, the A-terms in f l6.23ap - fl6.24cp can be eliminated by adding linear in z and z 
phases, which leads to 



U r . 



n e i4 



where 



(f)(z,z) = f Ea\ a (P,Q) z + / _1 eBAs(P,Q) z + constant. 

Using the definitions of W n and T n , and introducing h = h(Q)/h(P), 
f Q f p 
h = -^%, 

F P F Q 

we can write the final expression for the quasiperiodic solutions as 



(z, z) = u^ z ^ h 



eU n (z,z) + r A 

n \ 



(6.30) 
(6.31) 

(6.32) 
(6.33) 



e(c»(M)+/lv 

This time, the link between u n and r n is linear in terms of C, n : the transformation 
r n u n is achieved by C n ->■ C + f P UJ - 

The simplest of the quasiperiodic solutions f )6.33p . after some slight modifications, 
can be rewritten cnoidal wave: 

u n {z, z) = u,e iMz ' 2) sn ( n (z, z) (6.34) 
where the phase <p n and the function ( n are given by 

4> n (z, z) = £ z + r] z + 5 n, (6.35a) 

Cn(z, z) = i P z + 7] p z + 5 p n (6.35b) 



Explicit solutions for a (2+1)- dimensional Toda-like chain. 



20 



and snz = sn(z, k) is the elliptic sine. The parameters £ , p , Vo,p, ^o,p an d the elliptic 
modulus k are related by 

£ p r)p = sn 2 5 P (6.36a) 

ZoVp + £pVo = 2 sin 5 sn 8 P (6.366) 

and 

Corjo = \du5 p - e iSo cn5 p \ 2 (6.36c) 
for real 5 and S p . 

7. Solutions of the two-dimensional Volterra equation. 

As was mentioned in section [2j the authors of [5J [6] derived general solutions of the 
2DVE in the case of finite chain. Here we would like to present several other classes of 
solutions for this system, namely ones that can be obtained from the results presented 
in this paper using (13.19d) and ( 13.1961) . 

7.1. Toeplitz solutions. 



As follows from (13. 19 aj) and f!3.196j) . the constants r* and f* as well as the phases z/Q 



(that we defined in (14. 9p . (I4.10p and (I4.35p . (I4.36[) ) disappear from the final formulae 
for a n and b n , that can be written as 

An+lpn- 1 A n ~ lR n 

° n AnRn ' n AnD 

in the infinite case (— oo < n < oo) and 



(7-i : 



_ M n+l p w-l , _ M n-l p n /y ry\ 

an ~ Aknk-1 ' °«- Akok-1 y'- A > 

"n^n rt n D n-l 

in the semi-infinite case (1 < n < oo). One can easily obtain similar solutions with C"- 
and (^-determinants which we do not write here. 

1.2. Soliton solutions. 

Making the shift A n — > AnH^H^ 1 and introducing the matrix B n , 

B n = A n Kp (7.3) 
one can presents soliton solutions for (I2.5P in the following form: 

(7.4a) 



det 1 1 + A„H 5 H^ 


det 


1 + B„| 


det |1 + A n H ? 


det 


1 - 





1 det|l+A n H„ 


det 


1 - 


1- B n Hj,| 


c* det 1 + AnH^Hy 


det 


1 + B n | 



(7.46) 

with c* = f/hp and matrices A n = A n (z,z) and defined in section In these 
expressions one can find a symmetry l/6 n = a n (£ — >rj) which is a manifestation of many 
symmetries inherent in the 2DVE. 
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7.3. Periodic solutions. 

It follows from Ij3.19a|) and fj319fr|) that the ph ases that one has to introduce to eliminate 
the A-terms in f l6.23ap -f l6.24cl) cancel themselves and formulae (16. 15 a[) . (16. 1 5 eh yield 
periodic solutions for the 2DVE. After the shift C n — > £ n — u one can write these 
solutions as 



g(c + />)fl(c + j>) 



(7.5a) 



b n = - -± / V - N (7.5b) 



0(C„)*(c„ + J?« 



where c* = f^nfe® with £ n (z, z) and / defined in section [61 One can see that of two 
points that parametrize solutions (16.331) . P and Q, only one is left and, again, one can 
find the symmetry linking a n and b n : b n oc a n (A+> B). 



8. Conclusion 



We have obtained three types of explicit solutions for the model (11. ip . To conclude, 
we would like to discuss a few questions that have not been studied comprehensively in 
this paper. The first one is related to whether these solutions are real or complex. All 
presented above are, in general, complex ones. From the viewpoint of applications there 
are situations when namely such solutions are needed. For example, in the model of the 
graphite-like magnetics with Heisenberg-type interaction discussed in [3] u n oc +io~n^ 
where are the components of er n , er n = (^Jn \ o~n \ &nj, which means that in general 
situation u n is complex. However, sometimes one may be interested in the real solutions, 
as e.g., in the theory of hydrodynamic-type systems [2]. So, one faces a natural question 
of restrictions that should be imposed on the parameters of solutions that ensure their 
reality. This is a particular case of a problem which frequently arises in the theory 
of integrable equations, especially in the case of (quasi)periodic solutions, and surely 
deserves special studies. 

Another range of questions is related to the symmetries of our problem. The 
symmetry 5 — > —5 in (I3.13P and the following equations, which is another form of 
n — > —n symmetry, is the simplest one. The symmetries that we saw in the previous 
section, l/b n = a n (£ — > w) and b n cx a n (A -h- B) are less trivial. At the level of the 
ALH, they describe the links between various Miwa's shifts. Since the ALH is out of 
the scope of this paper (though implicitly we used some of its properties), we do not 
discuss this topic further, noting only that studies in this direction may lead to some 
kind of nonlinear superposition formulae for our problem. 

Finally, we would like to repeat that from the viewpoint of the 2DVE or RTC our 
main result, solutions of (jl.ip . are nothing but solutions of the auxiliary linear problems 
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(I2.17p — (I2.19p . In other words, after some additional work (such as, e.g., extracting the 
the dependence on A) one can obtain from our results the Baker- Akhiezer function of 
the 2DVE or RTC, which is usually more complicated problem that to obtain solutions 
for an equation itself and which can be useful for applications. 



Appendix A. Toepllitz determinants. 

Applying the Jacobi identity 

AA^ = A£A£-A^A£ (A.l) 

where A is the determinant of a matrix, A 3 k is the determinants of the same matrix with 
jth row and kth column being excluded, etc, to F k + i(C) for (ji, j 2 ; k\, A^) being equal to 
(1, t + 1; 1, i + 1), (1, 2; 1, i + 1), (1, 3; 1, i + 1) and (1, £;1,£+ 1) one can obtain 

X? := At, F k i+1 (() - A* F}(C) + (A,*" 1 F* +1 (C) = (A.2 

Z k := kU FjgCC) - A k+1 F{(C) + (tf F* +1 (C) = (A.3 

$ ■= 00 d+AU - A k+1 d + F k {Q + (A k d + F k+1 (() = (A.4 

5 ■= F, fe +1 (C) - A* 0_F*(C) + CA}" 1 <9_F* +1 (C) = (A.5 
correspondingly. These identities, combined with 

A k := {A k ) 2 - AlX+i ~ Ar X A* +1 = (A.6 

lead to 

Wj := A*F* +1 - A k +1 F k - A^Fjg = (A.7 

Y k := AfFj+J - A k+1 F k +1 + (A k +1 F k+1 = 0. (A.8 
Expanding the above equations, 

W k = A k + (w k + ... + (-O^w* (A.9 

X{ =C^ k + ... + (-CY- 1 4 (A.10 

Y* =Cy k + ... + {-C,YA k (A.ll 

Z k =Cz k + ... + (-CY^z k , (A.12 



where 



with 



D+A^-Af-A^A^ (A.13 
D + A*" 1 • At x - A^A* (A. 14 

Ajf ■ At i - A{l 1 1 A{ +1 (A. 15 



w£ = - A £ fc <9_A £ fc + Ati<9_Ati + Aj _1 9_Aj +1 (A. 16 

x| = D_ Ati • A £ fc + At'Ati 1 (A. 17 

yj = D_ Aj • At 1 + AKA{ +1 (A.18 

z £ fc = L>_ At x • A^ +1 + A^At (A. 19 
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and introducing 

3 k :=d + Z k -l k = (A.20) 

3 k := <9_X| -I* = (A.21) 

one can derive by straightforward algebra the following identities: 

A£_! FJ|i' • A k + a£}F{) = A k 3 k - {d + A k e ) Z k - F k w k+1 (A.22) 

A^ 1 (D + F k +1 ■ A k + (A k e+1 F k ) = Afj}" 1 - (0+ Af) Z^ 1 - (FX (A.23) 

Aj_ 2 (D+ F* +1 • A k + A k+1 F k ) = Apt, - (d + A k ) Z{_i + Af^tx 

- F k _^ +1 + CFfcfz} (A.24) 

and 

Afe 1 (D- F k H ■ A k - A k +l F k+l ) = A k 3 k e +1 - (<9_A £ fc ) X k+1 - F\ +1 y k (A.25) 

At, (D_ F k +1 ■ A k - (AtlFt 1 ) = A k 3 k - (d_A k ) X k - (F^yt 1 (A.26) 
A k _ 2 (D_ F k ■ A k + At'Ft 1 ) = AfiU + AfZ k _ x - {d-A k ) X^ 

+ F}_!Sf - CFfc 1 ^ 1 ( A - 27 ) 

which means 

D + F{+f -Af + ^F} =0 (A28) 

D + F k +1 -A k + CA k +1 F k = (A29) 

£>+F* +1 • A k + A k+1 F k =0 (A.30) 

and 

D.Fjfi.Aj-Aj^F^ 1 =0 (A.31) 

D_ F} +1 • A k - (AtlFt 1 = (A.32) 

D_Fj • Aj + Aj- 1 F{ +1 =0. (A.33) 

These equations, when rewritten in terms of B k given by ( 14.30 . are nothing but (I4.6a0 - 
( I4.7c0 . while in terms of C k given by ( 14.170 they become 

= D + C k ■ A k + (Ct, A k +1 (A.34) 

= D + C k+1 ■ A k + A k tl Cti (A.35) 

= D + C*+ 1 ■ A k +1 + A^ 1 C{ (A.36) 

and 

= D_ C k ■ A k - (Atl Ctl (A.37) 

= D_ C k ■ A k +1 + Atl C k+1 (A.38) 

= D_ C k+1 ■ A k - A k +1 Ctl (A.39) 
The above formulae are enough to obtain the Toeplitz solutions presented in section HI 
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